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Abstract. Magnetization process of the mixed spin- 1/2 and spin-3/2 Ising-Heisenberg 
diamond chain is examined by combining three exact analytical techniques: Kambe projection 
method, decoration-iteration transformation and transfer-matrix method. Multiple frustration- 
induced plateaus in a magnetization process of this geometrically frustrated system are 
found provided that a relative ratio between the antiferromagnetic Heisenberg- and Ising-type 
interactions exceeds some particular value. By contrast, there is just a single magnetization 
plateau if the frustrating Heisenberg interaction is sufficiently small compared to the Ising one. 



1. Introduction 

The spin-1/2 quantum Heisenberg model on the diamond chain has enjoyed a considerable 
research interest since two unusual tetramer-dimer and dimer- monomer phases were theoretically 
predicted by Takano et al. [1] as a consequence of the mutual interplay between quantum 
fluctuations and geometric frustration. Another stimulus for investigating diamond chain models 
was triggered by a recent discovery of molecular solids B14CU3V2O14 [2], CU3 (003)2(011)2 [5], 
Cu3(Te03)2Br2 [I], which represent possible experimental realizations of the diamond chain 
compounds. Motivated by these experiments, several phenomenological models have been 
suggested and solved with the aim to bring insight into a frustrated magnetism of the diamond 
chain compounds (see [5 J and references therein). In the present work, we will develop an exact 
analytical treatment for the mixed spin-1/2 and spin-3/2 Ising-Heisenberg diamond chain that 
should clarify how the geometric frustration influences a possible magnetization scenario. 

2. Model and its exact solution 

Let us write the total Hamiltonian of the mixed spin-1/2 and spin-3/2 Ising-Heisenberg diamond 
chain as a sum over plaquette Hamiltonians 7i = J2k ^fc> where each plaquette Hamiltonian TLj, 
involves all the interaction terms included within fcth diamond-shaped plaquette 

H k = ^HSsfc-rSsfc+JitSffc-i+SffcX/l^ 

Here, ft z k and = (S%, 8%, S|) denote spatial components of the spin-1/2 and spin-3/2 operators, 
the parameter Jh > labels the isotropic antiferromagnetic interaction between the nearest- 
neighbouring Heisenberg spins and the parameter J\ > accounts for the antiferromagnetic 



Ising-type interaction between the Heisenberg spins and their nearest Ising neighbours. The last 
two terms determine the magnetostatic Zeeman's energy of the Ising and Heisenberg spins in an 
external magnetic field H\ and H^i oriented along the z-axis. 

Exact solution for the aforedescribed Ising-Heisenberg diamond chain can be obtained by 
combining three exact analytical techniques: the Kambe projection method [6], the generalized 
decoration-iteration transformation [7] and the transfer-matrix method [8]. First, let us take 
advantage of a commutation relation between different plaquette Hamiltonians [Hi,'Hj] = 
valid for each i ^ j to partially factorize the partition function 

JV 

2 = E T Hs k} exp(-PH) = ^11^ expi-pHk). (2) 
{fj-k} {fj-k} k=1 

Here, j3 = l/(k-QT) (k-Q is Boltzmann's constant, T is the absolute temperature), the symbols 
J2{n k } an d ^{Sfc} denote a summation over spin degrees of freedom of all Ising and Heisenberg 
spins, respectively. To proceed further with calculation, it is necessary to perform a trace over 
spin degrees of freedom of two Heisenberg spins from kih diamond plaquette as indicated by the 
symbol Tr^ shown on the rhs of Eq. @. This calculation can easily be accomplished by making 
use of the Kambe projection method [6], since the Hamiltonian 7ik can also be viewed as the 
Hamiltonian of the spin-3/2 Heisenberg dimer in the effective field H c r = H-^ — J\{^ k _ 2 + ^ k+l ) . 
Consequently, the complete set of eigenvalues corresponding to the plaquette Hamiltonian 
Tik can be expressed solely as a function of two quantum spin numbers 5tot and 5f ot , 

Ek(S to t, = -15J H /4 + JnS t ot{Stot + l)/2 - S* ot H eS - Hi(fi 3k _ 2 + /4fc+i)A (3) 

which determine the total quantum spin number of the spin-3/2 Heisenberg dimer and its 
projection towards the z-axis. According to the basic laws of quantum mechanics, the total 
quantum spin number of the spin-3/2 dimer gains four different values S^ot = 0, 1, 2, 3, while the 
quantum spin number 5 t z ot gains 2£tot + l different values 5 t z ot = — S to t, — SW+lj • • • > •S'tot for each 
allowed value of 5 to t- The energy eigenvalues (J3J) can be then used to obtain the relevant trace 
of the last expression on rhs of Eq. ([2]), whose explicit form immediately implies a possibility of 
performing the generalized decoration-iteration transformation [7] 

3 n 

Tr fc exp(-/3?4) = ex.p[(3Hi(fi^ k _ 2 + ^3fc+i )/ 2 l E E exp[-/3J H n(n + l)/2] cosh(/?ff eff m) 

n=0 m=—n 

= Aexp[(3Rf, z 3k _ 2f i z 3k+1 + m{vl k -i + Mffc+i)/2]. (4) 

Note that unimportant constant term exp(15/3Jn/4) has been for simplicity omitted from lhs of 
Eq. (H]) as it does not basically change physics of the studied system. It should be mentioned, 
however, that the decoration-iteration transformation must hold for arbitrary combination of 
spins states of two Ising spins an d ^3^+1 involved therein. This self-consistency condition 

unambiguously determines so far not specified mapping parameters A, R and H 

A = (W+W.W 2 ) 1 ^, PR = ln(W + W^) - 2 In W, (3H = (3Hi - In W + + I11W-, (5) 

which can be uniquely expressed through the functions W± = F(±l) and W = F(0), where 

3 n 

F(x) = E E exp[-/3J H n(n + l)/2] cosh[/3m(Jix + H n )]. (6) 

n=0 m=—n 

It is easy to verify that the transformation (jlj) establishes a precise mapping relationship between 
the mixed spin-1/2 and spin-3/2 Ising-Heisenberg diamond chain and its corresponding spin- 
1/2 Ising linear chain with the nearest-neighbour coupling R and the effective magnetic field H 



given by Eqs. ([5]) and ([6]). A direct substitution of the mapping transformation (jl]) into Eq. ([2]) 
actually connects partition functions of both these models through the relation 

Z((3,J H ,J h H ll ,H l ) = A N Z lsing ((3,R,H). (7) 

It is noteworthy that the partition function Z\ s \ ng of the spin- 1/2 Ising linear chain can simply 
be attained with the help of the transfer-matrix method (see for instance Ref. [Sj) and thus, 
our exact calculation is formally completed. In this respect, the mapping relation ([7]) can also 
be used to obtain exact results for the Helmholtz free energy T = —f3~ l \nZ and subsequently, 
both the sublattice magnetization toh = — ^ {jni^) anc ^ m i = ~~k (JlTf)' 

3. Results and discussion 

Now, let us step forward to a discussion of the most interesting numerical results obtained for 
magnetization process. Before doing so, it is useful to mention that magnetic properties of the 
frustrated Ising-Heisenberg diamond chain with other two particular values of the Heisenberg 
spins S = 1/2 and 1 have been detailed discussed in our previous work [5]. The main objective 
of the present work is therefore to investigate magnetization curves of the mixed spin- 1/2 and 
spin-3/2 Ising-Heisenberg diamond chain, which should serve for the sake of comparison with 
the recently published results [5] and would shed light on how the scenario of magnetization 
process depends on a magnitude of the Heisenberg spins. 

To simplify further discussion, let us first rescale all interaction parameters with respect to 
the Ising-type interaction J\ by introducing the following set of dimensionless parameters: a = 
Jh/Ji, h = Hi/ Ji = Hu/Ji and t = k^T/Ji as describing a strength of the geometric frustration, 
external field and temperature, respectively. The total magnetization [m = (mj + 2m-n)/3\ 
together with both sublattice magnetization mi and mn is plotted in Fig. [T] against the reduced 
magnetic field for several values of the frustration parameter a. As one can see, we have 
altogether found six different types of magnetization curves. If the frustration parameter is 
a < 1/3 [see Fig. [1(a)], there exists just a single intermediate magnetization plateau before the 
magnetization tends towards its saturation value. On the other hand, two or three intermediate 
magnetization plateaus can be observed provided that the frustration parameter is from the 
interval 1/3 < a < 1/2 [Fig. 1(b)] or a > 1/2 [see Figs. QJc)-(f)], respectively. It is quite 
obvious from this analysis that the geometric frustration can be regarded as the driving force 
behind appearance of the greather number of intermediate magnetization plateaus. 

It is worthwhile to remark that the identified magnetization plateaus appear at 1/7, 3/7 and 
5/7 of the saturation magnetization when normalizing the total magnetization with respect 
to its saturation value. All these fractional values are in accordance with the Oshikawa- 
Yamanaka-Affeck rule [9], which has been proposed as a necessary condition for the formation of 
quantized plateaus. As it could be readily understood from the relevant results for the sublattice 
magnetization, the spin arrangements emerging at intermediate plateau regions with the same 
value of the total magnetization might significantly differ. However, the precise nature of spin 
arrangements to emerge at different plateau regions requires a more comprehensive analysis of 
the ground-state phase diagram, which is left as a challenging task for our future work. 

It should be also mentioned that the stepwise magnetization curves with an abrupt change(s) 
of the magnetization observable at critical fields occur merely at zero temperature, while there 
are no real magnetization jumps at any finite temperatures. It turns out, however, that the 
magnetization curve exhibits a rather rapid continuous increase in the vicinity of critical fields 
whenever the temperature is set sufficiently small. Contrary to this, the magnetization curves are 
gradually smeared out when the temperature is raised from zero and consequently, the plateau 
structure gradually diminishes from the magnetization curves above t ~ 0.4. 

In conclusion, we have developed an exact analytical treatment that enables to study the 
magnetization process of the mixed spin-1/2 and spin-3/2 Ising-Heisenberg diamond chain. 
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Figure 1. The total magnetization m (solid lines) and the sublattice magnetization mi (dotted 
lines), me (dashed lines) as functions of the external magnetic field at very low temperature 
t = 0.01 for several values of the frustration parameter (a) a = 0.2, (b) a = 0.4, (c) a = 0.8, 
(d) a = 1.3, (a) a = 1.8, (b) a = 2.0. 



Within the framework of this exact method, we have found a rigorous evidence for multiple 
frustration-induced plateaus, whereas the number of intermediate plateaus basically depends on 
a strength of the frustration parameter a. It is noteworthy, moreover, that the developed exact 
analytical treatment can be rather easily generalized to account for the magnetization scenario 
of the mixed spin-1/2 and spin-S Ising-Heisenberg diamond chain, which will be dealt with in 
our subsequent more comprehensive study. 
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